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1 Introduction 



The idea of the description of 2d theories of conformal matter coupled to gravity through 
the simpler models of the lattice type, which fall into the same universality classes, sug- 
gested in 0, ^ led to an explosure of the interest to different matrix models in 1989, when 
it was understood that the most convenient choice of the lattice models is related to the 
triangulations of random surfaces (in spirit of Regge approach [0), which can be coded 
in terms of matrix models . More concretely, there was proposed a class of multi-matrix 
models whose (multi-)critical points correspond to minimal series of conformal matter 
coupled to gravity H, 0, ||, ||, |10], |ll], |T2|, The general description of the (double 
scaling) continuum limit of these models essentially involves the following two properties 
of the double scaled partition function: it should be the r-function of proper reduced KP 
hierarchy and it should satisfy the W^-algebra [0, |], |15|. The combination of these two 
properties leads to the string equation. 

Indeed, it was the string equation that was derived firstly, and the general proper- 
ties of the matrix models were guessed later from the properties of the string equation. 
Unfortunately, it is very difficult to prove them completely, as the the structure of the 
theory is very complicated when considering higher multi-critical points and/or multi- 
matrix models . From the other hand, let us point out that the determing properties of 
the multi-matrix model partition function are observed a/ier taking the continuum limit. 
Therefore, it is reasonable to find out new multi-matrix models with these properties at 
the discrete level. In fact, at the discrete level one usually has a property of manifest dis- 
tinguishing the integrability and the W^-invariance simultaneously. Due to this fact, one 
can effectively investigate the continuum limit by observing only ly-algebra (by modulo 
one subtlety, see below the subsect.5.2). As for the integrability, it is the standard thing 
that one can continue the model off critical point preserving this property. This contin- 
uation is rather natural and simultaneously gives rise to a proper regularization of the 
continuum theory (we also know this phenomenon in quantum integrable systems, where 
to regularize the theory one can put it on the lattice preserving the quantum integrability 
- see, for example, |TB[). 

All this was partially realized previously. Namely, there was proved that at the discrete 
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level the partition function of the Hermitian one-matrix model is the r-function of Toda 
chain hierarchy [|1^] and satisfies the Virasoro constraints [|T8|, |T9|. It was also shown |2y] 
that these constraints turn into the continuum constraints of the paper 0. Moreover, 
one can check [^] that in this continuum limit Toda chain hierarchy really turns into 
KdV hierarchy which corresponds to the double scaling limit of the Hermitian one-matrix 
model H H 

Unfortunately, the situation is not so clear in the multi-matrix model case. In the 
general description the standard multi-matrix models proposed in |]1T], |12|, |13[ have no 
any symmetry at the discrete level (though they still correspond to an integrable system, 
that is, to Toda lattice hierarchy [0, |2^, |2^). It does the proof of their l^-invariance in 
the continuum limit very difficult. There were also suggested other discrete multi-matrix 



models with the same integrability property and ly-invariance [^, |2^. Their continuum 
limit was investigated but the complete answer is not obtained up to now due to 
technical difficulties. 

At last, in the paper ^ a rather natural multi- matrix model was introduced, with the 
property of W^-invariance included by definition. This model has less trivial integrable 
structure and corresponds to multi-component KP hierarchy with a special reduction of 
AKNS type [^]. This model will be the main subject under consideration in the paper. 
As the corresponding partition function can be naturally presented in the form of the 
correlator in conformal theory, we call it "conformal multi-matrix model " (CMM). 

We consider the paper as one of the series of two papers devoted to the investigation 
of CMM. Another one p6| mostly discusses the the integrability, while the continuum 
limit is described in less details. The present paper thoroughly deals with the continuum 
limit of CMM. But, for the sake of completeness, we describe general properties of the 
discrete model. More concretely, in the section 2 the definition of CMM is introduced, 
and some of its general properties are discussed. Here we are interested only in the 
continuum limit of VT-algebra , the continuum limit of arbitrary multicomponent KP 
hierarhy of such type being out of the scope of this paper. In the sect. 3 the general 
approach to taking the continuum limit is discussed. It is considered for the simplest case 



of the Virasoro constraints in the one- matrix model in rather invariant terms lEDH, and 



some natural generalizations are proposed. It turns out that this invariant formulation 
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naturally continues to the general case due to the very special form of the discrete W- 
algebra resembling the structure of the continuum ly-algebra of the paper (indeed, just 
the "conformal" form of the partition function in this model is in charge of this). In the 
sect. 4 the simplest non-trivial case of two-matrix model (ly'i^^^-constraints) is considered 
in full details, and the general case is discussed in the sect. 5, some important but tedious 
calculations being shifted to the Appendix. The conclusion contains the list of the main 
results as well as some general discussion. 

2 Conformal mult i- matrix models 
2.1 Formulation of the model 

To begin with, we would like to introduce conformal multi-matrix models in accordance 
with the paper [|I| (see also |^). First, we show that the simplest example of discrete 
Hermitian one-matrix model can be easily reformulated in these terms. 

Indeed, Hermitian one-matrix model {p = 2) can be defined as a solution to discrete 
Virasoro constraints: 

oo 

Ln = Y.ktkd/dtk+n+ E d^dtadU (1) 

fc=0 a+b=n 
dZ2,N/dto = -NZ2,N 

The Virasoro generators (|l|) have the well-known form of the Virasoro operators in the 
theory of one free scalar field. If we look for such solution in terms of holomorphic 
components of the scalar field 

[Jn,Jm] =n6n+m,0, [g,P] = l (2) 

the procedure is as follows. Define vacuum states 

JfciO) =0, (iV|J„fc = 0, k>0 

p\0) =0, {N\p = N{N\, (3) 
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the stress-tensor 

T{z) = l-.mz)]': = J2T^z—\ T„ = 1 5: :J_,J,+.: + ^ ^ J^Jb, (4) 

^ k>0 ^ 

a,6>0 

Tn\0) = 0, n > -1 (5) 

and the Hamiltonian 

Hit) =^T.tkJk=i V{z)j{z) 

y(^) =E^fc^', 3{z) = ^d<t>{z). (6) 

fc>0 V/ 

Now one can easily construct a "conformal field theory" solution to ([|) in two steps. First, 

L„(iV|e^W... = (iV|e^WT„... (7) 

can be checked explicitly. As an immediate consequence, any correlator of the form 

(iVle^WClO) (8) 

(iV counts the number of zero modes of G) gives a solution to (|1|) provided 

[Tn, 6-1 = 0, n>-l. (9) 

Second, the conformal solution to (|) (and therefore to ([|)) comes from the properties of 
2d conformal algebra. Indeed, any solution to 

[T(^),G]=0 (10) 

is a solution to (^, and it is well-known that the solution to (|10D is a function of screening 
charges 

Q^= fj^= fe^^'t (11) 
With a selection rule on zero mode it gives 

G = exp Q+ j^Q^ (12) 

Of course, the general case must be G ~ Q^^^^ ■ Nevertheless, this choice of one of the 
two possible screening operators has a clear algebraic sense which we will discuss below 



in this section. It can be justified by the special prescription for integration contours, 
proposed in [Q], which imphes that the dependence of M can be irrelevant and one can 
just put M = 0). In this case the solution 

^2,iv[t] = (iV|e^WexpQ+|0) (13) 
after computation of the free theory correlator gives well-known result 
Z2,N = {N^V / n exp (- tkzl) A%{z) = 

i=l 

= {NNo\ U{N))-^ J DM exp (- J2 hM^) (14) 

N 

=\{{Zi- Zj) 

i<j 

in the form of multiple integral over spectral parameters or integration over Hermitian 
matrices. 

In the case of p = 2 (Virasoro) constraints this is just a useful reformulation of the 
Hermitian 1-matrix model. However, in what follows we are going to use this point of 
view as a constructive one. Indeed, instead of considering a special direct multi-matrix 
generalization of (|T^ |TI|, |T2|, |I^] one can use powerful tools of conformal theories, where 



it is well known how to generalize almost all the steps of above construction: first, instead 
of looking for a solution to Virasoro constraints one can impose extended Virasoro or W- 
constraints on the partition function. In such case one would get Hamiltonians in terms 
of multi-scala.T field theory, and the second step is generalized directly using screening 
charges for W^-algebras. The general scheme looks as follows 

(i) Consider Hamiltonian as a linear combination of the Cartan currents of a level one 
Kac-Moody algebra Q 

i7(t«,...,t(-°'^^))= 4^Va^., (15) 

A,fc>0 

where {/Xj} are basis vectors in Cartan hyperplane, which, say for SL{p) case are chosen 
to satisfy 

1 P 
P j=i 

(ii) The action of differential operators VT/"'' with respect to times {t\!'^} can be now 
defined from the relation 

wf)(Ar|e^«*»... = (Ar|e^«*»wf^.. , a = 2,...,j9; i > 1 - a, (16) 



where 

W('^)(;^) =Y.--[f^xdH^)r: + ... (17) 

A 

are spin-a W-generators of Wp-algebra written in terms of rank ^-component scalar fields 

(iii) The conformal solution to ly-constraints arises in the form 

Z^Tlit}] = (Ar|e^«*»G{g(-)}|0) (18) 
where G is an exponential function of screenings of level one Kac- Moody algebra 

Q(") = ^jW = ^e"^ (19) 
{ck} being roots of finite-dimensional simply laced Lie algebra Q. (For the case of non- 



simply laced case see [28|. Below Q = SL{p) if not stated otherwise.) The correlator ([T 



is still a free-field correlator and the computation gives it again in a multiple integral form 

z^^m] - / n 



i=l \ \,k>0 

X n lim^'-^-^fr^ m 

{a,l3) i=lj=l 

The expression ( [20| ) is what we shall study in this paper: namely the solution to discrete 
VT-constraints which can be written as multiple integral over spectral parameters {-z^"^} 
(this integral is sometimes called "eigenvalue model"). The difference with the one-matrix 
case (|T^) is that the expressions (^) have rather complicated representation in terms 
of multi-matrix integrals. Namely, the only non-trivial (Van-der-Monde) factor can be 
rewritten in the (invariant) matrix form: 

f (a) 



n 11(4"^ - zfY^ = fdet{M(-) Af^'^)}]"^ , (21) 

j=ij=i 



where / is the unit matrix. Still this is a model with a chain of matrices and with closest 
neibour interactions only (in the case of SL{p)). 

Now we would like to say some words on the general structure of this model. Let 
us point out that its partition function is nothing but the correlation function of objects 



which have clear algebraic meaning. Indeed, the time dependent exponential is generated 
by Cartan currents of SL{p), and the exponentials of the screening charges correspond to 
the exponentials of other (non-Cartan) generators of SL{p). But we should stress that 
again it should not be an arbitrary combination of these generators, but only of any p — 1 
of them. Indeed, the different choices of these p — 1 generators correspond to the different 
models but with the same properties of the integrability and the l^-invariance . 

In fact, it is very immediate thing to fermionize these expressions and merely write 
down the expression for proper r-function of reduced KP hierarchy of the generalized 
AKNS type P^ , 3U| in the form of fermionic correlator This generalized AKNS 

reduction has its origin in the reducing from the general case of GL{p) algebra to its 
simple subalgebra SL{p) [^, when considering proper ly*^^^- algebra p7[]"^. 

The purpose of this paper is to show that CMM defined by (|2|) ClS db solution to 
the ly-constraints possesses a natural continuum limit. To pay for these advantages one 
should accept a slightly less elegant matrix integral with the entries like 



2.2 On the proper basis for CMM 

Now we would like to discuss briefly the manifest expressions for constraint algebras in 
terms of time variables. 

The first non-trivial example (which we use as a demonstrating example in the section 
4) is the p = 3 associated with Zamolodchikov's VFs-algebra and serves as alternative 
to 2-matrix model. In this particular case one obtains 

H{t,i) = 4= T.(tkJk + hJk) (22) 
v2 fc>o 

oo 

Ln = ^{ktkd/dtk+n + kikd/dik+n) + 

fc=0 

+ (dVdt^dh + dydtadh) (23) 

a+b=n 

= (ktkltid/dtk+n+i - kikltid / dtk+n+i - 2ktkliid / dik+n+i) + 

k,l>0 

Indeed, all this has an interpretation immediately in terms of proper Hamiltonian reduction - see also 




8 



' ' n 



+2E 



fc>0 



a+b=n+k 



+ 



+\ E {d''/dtadm,-d''/dtadhdU), (24) 

a+b+c=n 

where times tk and ik correspond to the two orthogonal directions in S'L(3) Cartan plane. 
(We use the standard specification of the Cartan basis: e = cti/v^, e = \/?)U2/ In 
this case one has six screening charges Q*-^"*'' {i = 1,2,3) which commute with 

W^'\z)=T{z) = ^:[dcf>{z)r: (25) 

and 

W^'\z) = j2:{f,,dcf>{z))': , (26) 

A=l 

where fi^ are vectors of one of the fundamental representations (3 or 3) of SL{3). 

This basis was originally used in |T]] as it just corresponds to the integrable flows, for 
the continuum limit we will use another basis in the Cartan plane connected with t ±t. 
(In other words, this is the question what is the proper reduction, or what combinations 
of the "integrable" times should be eliminated.) 

To begin with, we consider the simplest non-trivial case of p = 3. Then introducing 
the scalar flelds 



k k C'tfc 

d<p^'\z)=j:kt^^^z'-'+j:A)'-'-\ (28) 

k k Ot^. 

with t^^^ = [iik + tfc)/2v^, t'^^^ = [itk — tfc)/2v^, one obtains the expressions: 

W^^\z) = -:d<j)^^\z)d<j)^^\z): , (29) 

W^'\z) = ^Y.-\d<P^''i^)f-- ■ (30) 

instead of (H) and (|2|). 

This choice of basis in the Cartan plane is adequate to the continuum limit of the 
system under consideration, as the latter one is described by completely analogous expres- 
sions 0. Now let us describe this basis in more invariant terms and flnd the generalization 
to arbitrary p. 
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Comparing (|30|) with (p6[) , we can conclude that 90^*^ = /3id4> corresponds to the basis 

/3l,2 = ^(V^^2±^«2). (31) 

This basis has the properties 

f3,-(3, = l, f3,-(3, = 0, f3,-f3, = 0. (32) 

Now it is rather evident how this basis should look in the case of general p. Due to |0] 
we can guess what is the choice of the proper scalar fields: 

a0«(.) = E + E -^z-'^-\ (33) 

k k Otj^ 

This choice certainly corresponds to the basis with defining property (it can be observed 
immediately from the relations (|53|) and 



/3, • /3,. = (34) 



the proper choice of the Hamiltonians in (15) being 



H = Y.t^f3rJk, (35) 

i,k 

what determines new times adequate to the continuum limit. 

Let us construct the basis ( ^4]) in a manifest way. To begin with, we define a set of 
vectors {/Xj} with the property: 

^ i 



The iy(")-algebra can be written in this basis as follows ||27l 



l^(") = (-)"+i E n = 2,...,p. (37) 

l<ii<...<ife<p m=l 



Now the basis ( |5^ can be constructed from ( |5BD by diagonalization of the following cyclic 
permutation p|, |33| : 

Hi^tii+i, Mp^Mi i = l,...,p-l. (38) 

This transformation has as its eigenvectors, their manifest expressions being of the 
form: 

1 ^ 27ri 

/3fc = — Eexp{ jk}^ij, k=l,2,...,p-l. (39) 

10 



It is trivial to check that the properties (|34[) are indeed satisfied. One can immediately 
rewrite the corresponding VT-generators in the basis of /3j's. After all, one obtain the 
expressions similar to the continuum W^-generators 0, but with the scalar fields 
defined as in (|33D and without the "anomaly" corrections appearing in the continuum case 
due to the twisted boundary conditions. These corrections can be correctly reproduced by 
taking the p-th root of the partition function as well as simultaneously doing the reduction 
(see the sects. 4-5). 

Thus, the proposed procedure allows one to take the continuum limit immediately 
transforming the scalar fields as elementary building blocks. Neveretheless, before imme- 
diate doing of the continuum limit for any CMM let us consider this in the simplest case 



of the Hermitian one-matrix model EO] and get some insight for the general case. 



3 Double-scaling limit of CMM: preliminary com- 
ments 

3.1 Results of [20] for the one- matrix model 



To begin with let us briefiy remind the main points of pO|. 

It has been suggested in that the square root of the partition function of the 
continuum limit of one-matrix model is subjected to the Virasoro constraints 



£cont^/^ = 0, n> -1, (40) 

where 

E ^ + ^ Wt^;^^ + G E 

u-n \ ^/ Ol2(k+n) + l n^u^„ 



r 



cont 

fc=0 



dT2(k+n)+l 0<fe<n-l 

(n-fc-l)+l 



16 (16G) ^ ' 



are modes of the stress tensor 



where 



9^ = E ((^ + ^) 7^2„H-i^"-^ + ^^^"""^) • (43) 
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It was shown in |20| that these equations which reflect the VT^^-'-invariance of the 
partition function of the continuum model can be deduced from analogous constraints in 
Hermitian one-matrix model by taking the double-scaling continuum limit. The procedure 
(generalized below to CMM) is as follows. 

The partition function of Hermitian one-matrix model can be written in the form 



Z{h} = jvM exp Tr ^ t^M^ (44) 

•' fc=0 

and satisfies []TH|, ^ the discrete Virasoro constraints (|I|) 

L^Z = 0, n > 

fc=Q ULk+n 0<k<n '-"'k'^'^n-k 

In order to obtain the above-mentioned relation between VT-invariance of the discrete 
and continuum models one has to consider a reduction of model (^4D to the pure even 
potential t2k+i = 0. 

Let us denote by the r^*^ the partition function of the reduced matrix model 

r^^'{t2fc} = / I?Mexp Tr J] hkM^" (46) 

k=0 

and consider the following change of the time variables 

/ xn-mp (n+^) a-"-5 
n>m {n - my.T [m + 

where Qm = Tnt2m and this expression can be used also for the zero discrete time qq = N 
that plays the role of the dimension of matrices in the one-matrix model. Derivatives with 
respect to transform as 

« -g /C'tir* (48) 



dhk t'o{k-n-l)\T(n + l) df^n+i ' 



where the auxiliary continuum times T2n+i are connected with "true" Kazakov continuum 
times T^n+i via 



T2k+l — T2k+1 + , 1 /o ^2(fc-l) + l; (49) 



k 

'A; + 1/2 

and coincide with T2n+i in the double-scaling limit when a 
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Let us rescale the partition function of the reduced one-matrix model by exponent of 
quadratic form of the auxihary times T2n+i 

1 



r = exp -- AmnT2m+lT2n+l T\ 



red 

N 



m,n>0 



with 



A. 



r(n + |)r(m + |) 



\n+m^—n—m—l 



2P i 



n\m\{n + m + l){n + m + 2) 
Then a direct though tedious calculation demonstrates that the relation 



n+1 rrcd^red 
i+l-p -^2p ^ 

-red ' 

p=0 



Ec^+i(-ir 



is valid, where 



fe=0 



d 



+ E 



92 



dt2{k+n) Q<k<n 9t2kdt2{n-k) 



and 



k>l 



d 



-2k+l 



k>0 



^0 = E ( ^ + 2 ) 



C^^2(fc-1)+1 
9 



16^ 



fc>0 



^n= E ( ^ + 2 ) 



dT2k+i 
d 



+ E 



C^^2{fc+n)+l 

d d 



Here 



0<fc<n-l dT2k+l dT2{n-k-l)+l 

are binomial coefficients. 



> 1. 



(50) 



(51) 



(52) 



(53) 



(54) 



" p\(n-p)\ 

These Virasoro generators differ from the Virasoro generators (|4T| ) ^ by terms 
which are singular in the limit a — > 0. At the same time U^p'^^'^'^ the r.h.s. of (|52D do 
not need to vanish, since 

<2fc+i=0 ' ~~ <J''^2i+ldt2{n-i-l)+l i2fc+i=0 



= L2pT 



(55) 



^ j dt2i+ldt2{n-i-l)+l 

It was shown in that these two origins of difference between (5T) and ( |5^ actually 
cancel each other, provided eg. (^2]) is rewritten in terms of the square root \/f rather 
than f itself: 



n+1 

'"Ec^+i( 

p=0 



r 



*2fc+l— 



;i + 0(a)). 



(56) 
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3.2 Generalization of the Kazakov variables 

Now we would like to generalize the procedure proposed in the sect. 3.1 to the case of 
the general multi-matrix model . However, before doing this, we shall describe a simpler 
way to guess proper Kazakov variables as well as tilded time variables and the connection 
between discrete and continuum VT-algebra . It can be done already at the level of the 
"leading" (or " quasiclassical" in accordance with |jl8|) terms, i.e. those including only the 
first derivatives. It is possible to do due to the fact that the constraint algebra of CMM 
strongly resembles the structure of the continuum W-algebra [2] (see (^91) , (|30D ). 

To begin with, we would like to consider the simplest multi-matrix case of two matrices 
in details. The general case which is completely analogous to the two-matrix one but 
requires considerably more tedious calulations will be considered in the Appendix. 

The leading terms of the W^^^ and W*-^-' generators are given by the formulas 

d_ 

i=l k,l>0 (^4+n+« 



^^(3)lead ^ ^ ^ kt^it?-^, = ktP, (57) 

Oil 



d 

' 3m+l" 



+ + T3fc+2T3^+2-~ . (58) 

Now, to connect these two algebras, we have to find out a change of time variables similar 
to (0), (^Sp which gives the relation between leading terms of the generators (p7[), (pB]). 
This problem can be solved in two steps. First, we make a reduction of the discrete times 

4ti=4t2 = 0, ^ = 1,2, k = 0,l,... . (59) 

Second, we note that the changing of time variables 

Tfn + l + i^a-^-t 

:1m ^ 



9ln^ = E : ^3n+, , ^ = 1, 2, (60) 

n>m [n - my.T [m + 



- E T, ^—zS ^ = 1' 2' (61) 
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where is equal to mt^^l^ 0, and the auxihary continuum times T^n+i are connected with 
the Kazakov continuum times T^n+i through 



k 

'A; +7/3' 



Tsk+i — T^k+i + I „■ /o ^3(fc-i)+i- (62) 



It gives rise to the following relation between w'^^^'"^^'^ and y^'^^^'"^^'^ 

n+2 

E C;+'(-)"~Wif = >Vf 1^"^ n > -2, (63) 

p=0 

where VV^^)'^^^ is given by ( ^8]) with shifted limits of the summation 



fc+m>-n ^ ^ C'J 3(fc+m+n)+2 

+ E (/c + l) (m+|)T3fc+2T3^+2-~ ^ . (64) 



The summation in (plj) includes the terms with times T^, m < defined by trivial 
continuing the equation (pO|). 

The proof of (B3) is based on the formulas (^), ( pT| ) and the following identity for the 



F-functions 

i-Vn^ + h) ^ T{h)T{c-h + a) 
^o7!(«-7)!r(c + 7) r(a + l)r(c + a)r(c-6)- ^ > 

To eliminate the incorrect negative mode terms in ( p^ we use the partition function 

rescaled by an exponential of a quadratic form of the auxiliarly times T^n+i (compare with 

m) 

f = exp - Yl ^mnT3m+lT3n+2 Tn^,N2^ (66) 
\ m,n>0 j 

where is the partition function of the conformal two-matrix model after the reduc- 

tion (|59| ) and 

^ ^ r (n + I) r (m + I) 

r(i)r(|) n!m!(n + m + l)(n + m + 2)' ^ ^ 



^Similar to the case considered in the previous subsection we have to include properly into the set of 
the independent discrete variables the sizes of matrices Ni. These are just new variables of zero indices 
g'o' . In the continuum limit they are related to the cosmological constant (generally there should be 
many different cosmological constants). All this can be treated completely along the line of ||2^, so we 
will not discuss it here. 
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This rescaling of the partition function is equivalent to the following transformation 
of the generators VV^'^^ 



y m,7i>0 J \m,n>0 

where now not only leading terms should be taken into account in (|68|). 



Then direct calculations similar to the calculations in ||20[ show that the changing of 
the times (|60|), ( |6l]) and the rescaling (^) establish the relation like (|63|) 



T 



p+2 
n=0 



\p~n on 



(3) red red 



r 



red 



-, P>-2, 



(69) 



where the generators W^^ ^^'^ depend only on the times t^''^^' and derivatives with respect 
to them, and the generators VV^^-* are exactly the generators of the Vr*^^''-symmetry of the 
continuum model 0. 

By definition, the generators VF^^-* annihilate the partition function (PO) for p = 3 
while W^"^ ^'^'d^red ^Qgg j-^Q^ need to vanish, since 



(1,2) 



=0; i=l,2 



(3) red red 



3n 



r'=" + . . . 



(70) 



where . . . means terms similar to J2m.,k 9^ logT/<9t3m+iC^^3n-3(m+fc)+if^^3fc-2 with all possible 
correct gradations which do not vanish under the reduction (|59|). In analogy with the 
results of [20|, these additional terms provides the eq. (|69|) can be rewritten in the terms 
of the cubic root of the r-function (see the sect. 5. 2 for detailed explanations) 



yy{3) cont 3/^ 



p=0 



-p ^'^Sp ' 



T 



(l + 0(a)). 



(71) 



We exclude here all simple but tedious calculations of general case, since for our 
consideration of CMM below we will use a more efficient way to deal with the changing 
of the time- variables t ^ T (which was also proposed in |2^), that is, a scalar field 
formalism. For example, in the case of the one-matrix model the Kazakov change of the 
time variables (|47|) , (^8]) can be deduced from the following prescription. Let us consider 
the free scalar field with periodic boundary conditions ((|33|) for p = 2) 



difiu) = J2 9kU 

k>0 



2k-l 



E 



d 



k>l ^'^'ik 



U 



-2fc-l 



(72) 
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and analogous scalar field with antiperiodic boundary conditions 

d 



z ^ i 



2k+l 



Then the equation 



\d^z)f = a^dip{u)T 



red 



u = 1 + az 



(73) 



(74) 



r r 

generates the correct transformation rules (^Tf), (^) and gives rise to the expression (^ 
for Anm- Taking the square of the both sides of the identity (0), 



lT{z)f = -^T(M)r 



red 



(75) 



one can obtain after simple calculations that the relation (|52D is valid. 



4 The case of VF^^) in scalar field formalism 
4.1 Scalar field formalism 

In this section we would like to consider the case of iy*^^)-algebra in full details in the 
framework of the scalar field formalism described at the end of the previous section. In 
this subsection we describe the formalism, and, in the following ones, we apply it to 
Virasoro and iy*^^)-algebra respectively. 

Thus, let us consider the set of scalar fields for the discrete two-matrix model 

d 



fc>0 k>l Oty.k 

d_ 



dv''\n) = E + E 771) (76) 

A:>0 k>l Ut^^ 



and the scalar fields of the continuum model 

d 



d^^^){z) = Y.((k + \)T,,^,z^-i + 

k>0\^ ^J-3k+2 



-Z ^ 3 



9$(^)(.) = E((a: + ^)t, 



\k-,2z'-'^ + -^z-'^-U . (77) 



Then all the relations (pDj), ( [HT| ) and ( p7D can be encoded in the equations 

i^^W {z)f = au'-^^dif'^^ {uy'^, u^ = l + az, i = 1, 2. (7J 



J- • • ^red 
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Let us prove (^) for the case of i = 2 only. The i = 1 case can be treated analogously. 
We have to check that the relation = 1 + az, with using ( |60D and (|6TD, gives rise to the 
following equation 



k>0 k>l Ot^k 



= E + I) Tsk-,2z'~^ + a E ^^-'-'^ - E Amkfsn.^2Z-'-^ . (79) 

fc>0 ^ k>0 yO'-'sfe+l m>0 J 

Let us start our check with the derivative terms. Using the obvious relations 



u'- = (1 + azr = E ^ ("^)"""^ (80) 

™>o (a + 1 - m) 



and (|6TD , we obtain 



E ^ (1) 

fc>l Ctg;, 



-3fe-l 



^ g ^ r (fc + 1) r (-A: + 1) af—'^ 



fe>i /=o m>0 



(A; - / - l)!r (/ + I) m\T [-k - m + f ) ^Tsz+i 



af d 7rsinf7rA; + 4f ) a-™-'=-i , ^ 



/>o r (/ + I) 5^3/+i fci^o^'o (A: - /)!m!r [-k-m-\ 



'Vsin(7r7 + f) g ^ ^8l^ 

/t^o4^or(/ + |)r(-7-i)5%+/ /t^o^W'^ 

where we replaced the variable of the summation k ^ 'j = k + m and used the identity 

V ^ ; TT = ^/7- (82) 

The first sum in the l.h.s. of ([79|), after the substitution = 1 + az and using (|60D, 
falls into two different sums 

k r fn + I) (-)n-k^k-m-n-l 

l^^^zk^ - T^YT TTTT 3«+2^ 

fc>o fc>on>fcm=o r (^fc -m + Ij (n - fc)!m! 

k>on>km>k+i i (^/c — m+gjl^n — fcjlm! 
The first sum in (|83f) can be rewritten in the form 



r( 




1 ^ ^n— fe^fc— m— n— 1 


r 


{k~m + l) 


1 (n — A;)!m! 



E E E 'Trr, rr; — -^371+2^ 
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n Y (n -\- -] / 9\ 

EE^3n+2-) ^a^-"-V-^5,,, = a-i5: U + -)T3„+2^"-^ (84) 

n>07=0 r(7+o) n>0 ^ 



where we used the identity (|^). After introducing the new variable of the summation 
f3 = m — k — 1, changing the order of the summation and using the identity 

n (_\k 1 

^^{n-k)\{k + [3 + l)\ f3\n\{(3 + n+l)' ^ ' 

the second sum in the r.h.s. of ( PBD can be rewritten in the form 

r(^ + l) 



\n „—n—l—fi 



3 



-1 



^t^on'To ^'"^'r - i) (/?)!n!(/5 + n + + n + 2)- 



r(^ + i) 



'EE^-/3^3n+2^-^-3, (86) 
/3>0n>0 



and from the last two lines of (|86D the expresion ( |67D for the matrix Anm follows. 

Performing the similar calculations for the case of z = 1 in (|78D, one can find the same 
expression for the matrix Amn- So, we may conclude that the changing of times (|3)? (HD 
results in the equations (ffHI). 



4.2 Virasoro constraints of the two- matrix model 

It was proposed in that the continuum "two-matrix" model possesses the W*-^'' and 
Virasoro symmetries, the Virasoro generators £„ and the generators of the W^^^ algebra 
being constructed from the scalar fields ( |77D in the following way 



nz) = i:a$«(z)9$(^)(z): - ^ = E J?^, (87) 

W(^)(.) = ^ (: (9$«(.))': + : {d^^^^z))' -^J = E (§8) 

We will prove below that the equations ( |78D yield the relation between the generators 
Cn and the corresponding generators W^'^ ^'^'^ = associated with the reduction (^91): 

-1 n+l rred_red 

^Cr^f = a- E C;^\-r^'-'^^, n > -1, (89) 

p=0 

19 



where Virasoro generators Lg*^ are defined by the formula (|29|) , where only tg^ and cor- 



responding derivatives are nonzero and the generators 



to 



E 

fc>i 

E 

fc>0 



d 



3fc+l 



- T.: 



d 



+ {k + 



d 



3fc+2 



3fc+l 



d ^ 



-T1T2 



' 3fc+2 



3k+2, 



E 



k—m=—n 



+ E 



9 



3fc+l 



9 



C^3m+1 

d (-)" 



9 



3 

n > 1 



' 3fc+2 



ST.: 



3m+2 , 



(90) 



a+fc=n-l C^^3m+1 9a" 

differ from the continuum generators £„ of the ref. by singular c-number terms. Indeed, 
again Lf^^s do not annihilate t^"^"^ exactly, but these two effects cancel each other, provided 
eq. ( p9D for f is rewritten in terms of the cubic root of the r-function In other words, 
doing accurately the reduction procedure in the Virasoro constraints of the discrete two- 
matrix model one can rewrite (^) in the form 



1 fo T 



n > -1. 



(91) 



t3k+i=0; i=l,2 

Thus, we conclude that, indeed, the continuum Virasoro constraints for the case of p = 3 
can be derived from the corresponding Virasoro constraints of the discrete conformal 
two-matrix model. 

To prove (^) let us calculate 

<l>«(2)f $(2)(2)f 



E- 

ra>0 



-n-2 



91nr 



91nr 



' 3fc+2 



A:>0 

E 



91nr 



dlnf 



+ 



+n)+2 k=0 (^3k+2 dT3(^n~l-k)+2 \ 



+ - 



k>l 



91nr 



+ k + 



dlnf 2 

-'3A:+2 — ~ ~)~ — -'l-'2 

5^3(fc-l)+2 9 



^^3(^-1)+! 

where in (^2]) only singular at z —>■ oo terms are taken into account. 
Using at n > 1 the identity 

r(fc + |)r(m+|) (-)^-ia-"_ 



(92) 



\n-l 



E A 



km 



E 



fc + m— n— 1 
k,m>0 



fe+m — n— 1 
A;,m>0 



fc!m!r(|)r(i) 



n(n + 11 



9a" 



(93) 
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the formula (p^ ) can be rewritten in the compact form 



n>-l 



-n-2 



d'^ Inr 



r 



3m+l 



(94) 



Using the fact that the generators U^^ of the discrete Virasoro algebra can be obtained 



from the formula 



(95) 



one can easily show that 



u 



T 



n>0 



E 



d"^ Inr 



;)^(2).(i) 

m,fc>l 



(96) 



where again only the terms singular at m ^ oo are taken into account. Now using the 
equations (p8[) we can conclude that 

d'^ Inr 



E 

n>-l 



n-2 



T 



- E 



k+m= 



dTsk+2dT, 



n>0 



-3{n+l) 



-^3nT 



3m+l 



9^ Inr 

^ ;.,(2) (1) 

m.fc>l 



E 



(97) 



It follows now from (|80|) that 



jC„t 



E 



9^ Inr 



fc+m — n — 1 f^3fc+2f^r3m+l 

fe,m>0 

If we take into account that 

n+l 

a-^(_i)"+i-pc;+^ E 

p=0 



n+l 

a-"E^n+i(-)"+'-^ 

p=0 



d"^ Inr 



n.{2),(l) 



E 



(98) 



92 



m, fc> 1 



n+l p-1 

xEE(-irCn-fi 



/?oS ^3/+i5f3,+2 r (/ + 1) r + 
r (^m + 1) r (^p - m + 



p=2m=i (m-/- l)!(p-m-^- 1)! 

we may conclude that (|89|) is correct. 



E 



92 



f,9>0 



(99) 
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4.3 Connection between and Vl^? 



Now we are going to repeat the procedure of the previous subsection for the ly^^^^-algebra 
generators. The connection between the discrete and the continuum cases was formulated 
in the section 4.1 and based on the formula (^) which can be easily proved by means of 
the equations (|78D . Indeed, it is easy to see that the generators Wj^^^ '^^'^ of the reduced 
discrete model can be rewritten using the scalar fields ([76|) as 



E 



u 



3n 



3^3 



(100) 



Then the relation between the generators of the W^^^ symmetry for the discrete model 
and the generators VV^^-' is 



1-r 



T = a 



T 



p+2 



p+2\ 



W(3) red ,ed 
\ p—n ' ' an ' 



n=0 



and follows from the identity 



r 



1 fdip^^\u)T 



red ^ 



^red 



p>-2 



I ^rcd 



(101) 



(102) 



The generators VV^^-* n > —2 are defined by the formulas (^) and the generators VV. 



(3) 



and W_i (but not the yV^') have additional terms, cubic in times: 



;(3) 



^Tl + ■■■ and W3 = + ^TlT, + 



(103) 



To obtain ( |101|) from (|102|) one need the set of the combinatorial identities 

n+2 ^3 ^3 



E 



E 



n+2 



+2l )" E ( fi'i+p-n+1 + „ ^z; , o,lz;o, 

/=0 rn,p>l \ Ot„_i_^_p j Ol^^Ot^p 

E W + i;)T3a+i— — + 



Slnr 



E 



2 <9T3^+29T3^+2 



n+2 

a ^ C;[_,_2( — )" ^ ^ I Qm+p-n 
f=Q m,p>l 



,(2) 



+ 



(91nr 



'-"'n-m-p I '-"'3m.'-"'3p 



E 

f3+-f—a=n 



-(1) 

+ E 



(1)55.(1) 

91nf 



2\ ^ 91nf d'^ 

a -\ — I T^a+2 — = ~ 

3/ (9T3/3+i9T37+i /3+7+a=„_l C^3a+1 ^^^3/3+1^^^37+1 
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which can be proved by the direct calculation using the formulas (|60D, (|6l|) and identities 
for the F-functions similar to 

Tim + b)T{c-m) ^ T jc - d)T {b)T jb + c) 
^0 m\{a-m)\ d\T{b-c-d) ' ^ ' 

5 General case 

5.1 The continuum limit of VF^"^^ -algebra 

It is easy to generalize the formalism of the scalar field described in the sect. 4. 2-4. 3 to the 
general conformal multi-matrix models using the scalar fields with Zp-twisted boundary 
conditions. Let us introduce p — 1 sets of the discrete times t^*"*, i = 1, 2, ... ,p — 1 and 
k = 0,1, . . . for the discrete {p — l)-matrix model and consider the reduction 

t^l+^ = 0, t,j = l,2,...,p-l, k = 0,l,... . (105) 

One can choose the "discrete" and the "continuum" scalar fields in the form 



k>0 k>l utpj. 



a4«(.-) = Y.\{k+'-\ r,.«.'-V + ^^.---^l , (107) 

fc>0 I \ "/ Cllpk+p-i ) 

Tpk+i = Tpk+i + a-^--j— y^Tp(fc_i)+j, z = 1, 2, . . . ,p — 1. (108) 



Then the equations 



au'-P+^-^df^'^ {uy^ = {z)f, uP = l + az, i = 1,2, . . . ,p - 1, (109) 



generate Kazakov-like change of the time variables 



■n- 

a p 



(-)"-™rfn + l+ ,^ 

9^^= E : \,^/ t = l,2,...,p-l (110) 

[n — m)lL im + -1 



n>m 



A^)-T.- . . . ^ . = l,2,...,p-l. (Ill) 



dt'^;> ^o{k-n-l)\T(n+l + ^] dTp^+r 
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Then the consideration similar to those used in the previous subsections shows that the 
relation between the tilded continuum generators yV^\ i = 2, . . . ,p of the W-symmetry 
and the reduced generators of the discrete VT-symmetry W^fj^'^'^ is as follows 

1 _ n+i-l TXA(i)red^red 

~ W»f = a- E C:^U-r^'~'-' , n>-^ + l, (112) 

where the rescaled r-function is defined 

p-i 



r = exp 



^ (■) ~ ~ Id 

2 51/ $Z ^mn^pm+iTpn+p-i j T {tpk} (-'--'-3) 



4=1 m,n>0 



and the matrices A^}^ are determined by 

A^l = — J ( ; ■^ " , ,r %7 ^, z = l,2,...,p-l. (114) 

r(^)r(2^) n!m!(n + m + l)(n + m + 2)' ' ' v ; 

The proof of the equivalence of ( |109| ) to (|110|) , (|111|) , ( p.l3| ) and (|114|) is the same as 
in the sect. 4. 2-4. 3 and we omit it here. To prove the equations ( p.l2| ), we need the exact 
expressions how the generators VV^*-' and wjf^'^'^'^ are connected with corresponding scalar 
fields (|106D, ([1071) . 

Performing the proper reduction procedure ( |105| ), which eliminates all time- variables 
excepting those of the form tpl (i.e. leaves the 1/p fraction of the entire quantity of 
variables) we can obtain the relation 



1 n+i-l W^^'^T 



*p;V0 only 



;i + 0(a)), n>-i + l, 

(115) 

where W^*-* 's are the W-generators of the paper . Thus, we proved the ly-invariance of 
the partition function of the continuum p — 1-matrix model and found the explicit con- 
nection of its partition function with the corresponding partition function of the discrete 
{p — l)-matrix model. 



5.2 On the reduction of the partition function 

To conclude this section we would like to discuss the problem of reduction ( |5UD and ( |1U5D of 
the partition function in detail, with accuracy up to (non-leading) c-number contributions 
and only after the continuum limit is taken. More precisely, we reformulate the condition 
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of a proper reduction in the continuum limit in order to reduce it to more explicit formulas 
which can be immediately checked. As a by-product of our consideration we obtain some 
restrictions on the integration contour in the partition function (pOl) . 

To get some insight, let us consider the simplest case of the Virasoro constrained 
Hermitian one- matrix model [^]. Before the reduction the Virasoro operators read as in 
(|45|) . Then their action on logr can be rewritten as 

d log r d log r 



^ d\ogT ^ d"^ log r 
After the reduction, we obtain 



+ E 



d log r 



red 



2k- 



E 



(9^ log r 



red 



+E 



dt2k+2n m 9t2mdt2n —2m 

d log r'''^'^ (9 log T^'^'^ 



E 



Q2 log^red 
dt2m+ldt 



0. 



(116) 



2n-2m-l 



+ 



9t 



2m 



9t 



2n—2m 







under the condition 



9 logr 



red 



odd 



0. 



;ii7) 



;ii8) 



todd— 



The last formula is a direct consequence of the "Schwinger-Dyson" equation induced 
by the transformation of the reflection M —>■ —M in (^If). Indeed, due to the invariance 
of the integration measure under this transformation one can conclude that the partition 
function (^) depends only on a quadratic form of odd times. 

Thus, the second derivatives of logr over odd times do not vanish, and are conjectured 
to satisfy the relation 

^2lQg^red ^ ^2lQg^red 



E 



E 



(119) 



dt2mdt2n-2m m C^^2m+lf^^2n-2m-l 

where the sign "~" implies that this relation should be correct only after taking the 
continuum limit. In this case one obtains the final result (cf. (|56D) 

(9 log \/t^ V- log 
l^kt2k^^ + 2^ 

k 



dt2k+2n m 9t2mdt2n-2m 



0. 



(120) 



Thus, it remains to check the correctness of the relation ( p.l9| ). To do this, one should 
use the manifest equations of integrable (Toda chain) hierarchy, and after direct but 
tedious calulations one obtains the result different from the relation (|119|) by c- 
number terms which are singular in the limit a ^ and just cancell corresponding items 
in ( |54D (this is certainly correct only after taking the continuum limit). 
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All this (rather rough) consideration can be easily generalized to the p-matrix model 
case. In this case one should try to use all W^^-'-constraints with 2 < i < p. Thus, the 
second derivatives should be replaced by higher order derivatives, and one obtain a series 
of equations like ( |116| ). It is the matter of trivial calculation to check that these equations 



really give rise to the proper constraints satisfied by ^ (cf. (^) and ( |115| )) provided by 
the two sets of the relations like ( |118|) and (|119|) . 

Namely, the analog of the relation ( |118| ) in the p-matrix model case is the cancellation 
of all derivatives with incorrect gradation, i.e. with the gradation non-equal to zero by 
modulo p. The other relation ( |119| ) should be replaced now by the conditions of the 
equality (in the continuum limit) of all possible terms with the same correct gradation. 
In the simplest case of p = 3 these are 

d'^logr'^'^ ^ 92 log r""^"^ 



— 3m,— 1 m —3m 
^Slog^red 93 log r"^^^ 
/ ~ / 



log T^'^'^ ^ log 



^ ^_^w6r ^ ^^21) 

m,k '^^3m,+2dt3n-3{m+k)+ldt3k ^ d't3mdt^n-3{m+k)d't3k 

Again, this second condition is correct modulo some singular in the limit of a ^ 
terms, which appear only in the case of even p. Unfortunately, we do not know the way to 
prove this statement without using the integrable equations, what is very hard to proceed 
in the case of higher p. 

On the other hand, the cancellation of derivatives with incorrect gradation can be 

trivially derived from the "Schwinger-Dyson" equations given rise by the transformations 

M exp |2^|m (0 < k < p) of the integration variable in the corresponding matrix 

integral, the integration measure being assumed to be invariant. In its turn, it implies that 

the integration contour, instead of real line, should be chosen as a set of rays beginning in 

the origin of the co-ordinate system with the angles between them being integer times — 

p 

|. This rather fancy choice of the integration contour is certainly necessary to preserve 
Zp-invariance of p-matrix model system. 



red 



■'Simultaneously it means that the size of all matrices in ( ^T| ) should be specially adjusted, say, to be 
times of p. 
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6 Conclusion 



To conlude, we have demonstrated that the partition function of a new class of multi- 
matrix models (CMM) proposed in can be described by the same constraint algebras as 
in the case of the standard multi-matrix models |Tl|, |1^, |13|. It was conjectured in [Q] that 
it defines the partition function in full. Therefore, one can conclude that the partition 
function of CMM lies in the same universality class (i.e. have the same (double scaling) 
continuum limit) as the standard models. Certainly, there is an essential difference in 
these models both in their formulation at the discrete level and in the procedure of taking 
the continuum limit. Indeed, in contrast to the standard multi-matrix models, CMM's 
satisfy l^-constraint algebra already at the discrete level. But as a price, it leads to very 
complicated (though still of the nearest neighbours) interaction of matrices in the matrix 
integral. From the other hand, the correct reduction procedure corresponds to taking 
the square root of the partition function in the standard case, in contrast to the root 
of the p-th degree, which should be involved for the proper reduction of the conformal 
{p + l)-matrix model. 

Now we would like to emphasize that the double scaling limit of the standard multi- 
matrix models was never investigated honestly in the whole space of possible coupling 
constants {i.e. not nearby a critical point). The problem certainly is that it is very 
difficult to describe any matrix model in the whole space without any explicit constraint 
algebra imposed on the partition function. Therefore, the continuum limit for multi- 
matrix models proposed in should be considered rather as a definition than as the 
theorem. 

In fact, we know another case when multi-matrix model possesses ly-invariance . It is 
the partition function described in |25], whose continuum limit was investigated in . 
The arguments mentioned above suggest that it should be possible to propose adequate 
procedure of taking the continuum limit for the VT-constraints in this well, but it 

is not done up to now. 

At last, we would like to stress some points which are to be understood better. 

1) In the paper we have proved the correspondence between discrete and continuum 
constraint algebras at the level of leading terms (see the sect. 3. 2), but the whole proof 
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taking into account all corrections is still unknown. It is not a problem to do this in any 
concrete case, and we have demonstrated it in the case of VT'^^^-constraint algebra, but it 
is very difficult to deal with the general case. One of the main reasons is that it is difficult 
to write down the general expression for the generators of W^"^ algebra (see, for example, 

ii)- 

Nevertheless, the main statement should be correct. We think that there should be 
more elegant and efficient way of doing the continuum limit, and that to find out this 
one should understand deeper the role which is played by the scalar fields with different 
boundary conditions. 

2) We do not know how it is possible to prove the condition of proper reduction like 
( |121| ). It was done explicitly only in the one- matrix model case (the relation (|119| )), where 
we used the equations of the integrable hierarchy (Toda chain hierarchy) in the manifest 
way and after rather tedious calculations proved the statement [^. But even in the 
case of two-matrix model the integrable equations are very complicated (maybe, we do 
not understand them completely) and it is absolutely unclear how one should solve this 
problem in the general case. Again, we think it might be done in a simpler way. 

The solution to all these problems seems to be important for a better understanding 
of the role of the double scaling continuum limit. We hope to return to all these problems 
elsewhere. 
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8 Appendix 

In this appendix we consider Kazakov-like change of the time variables in the general 
case of the conformal p — 1-matrix model as well as the relation between the discrete and 
the continuum VT-algebras which can be easily extracted from the analyzes of the leading 
terms only (see the sect. 3. 2). We can use again the fact that the constraint algebra of 
CMM, by construction, strongly resembles the structure of the corresponding continuum 
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constraint algebra. This assertion means that the generators W^^^^s of the discrete algebra, 
after the reduction ( |105| ), depend on the scalar fields ( |106| ) through 



sez " i=i 



(122) 



where the sign "~" means that we omit some numerical coefficient. Then it is obvious from 
(|122|) that the generators W^)'''s, are given by the sum (with proper binomial coefficients) 
of the terms 



d 



a=lfe„>0 6=1 mb>l (^tprrib 



1 otj 



V-3 3 
6=1 a=l 



(123) 



where j = 0, 1, . . . ,p — 1, i = 1, 2, . . . ,p — 1. 

Now we are going to prove that the Kazakov change of the time variables ( |110| ), (|111 
results in the relation 

p-j 



n+p—l 
s=0 



a=lka>0 b=l mi,>l Otpmi, 



mf,=Y^ ka+i~p+s+l 



nn E T.[na+-]T,^. 



d 



a=lfe=l /i,>0 

y , na>-n-i+p-j 



Pfb+p-i 



(124) 



J2 ■fb=J2 na+n+i-p+j 

where the sums over a and b in ( |124D run in the same limits as in ( p.23| ) and again in the 
sums over some incorrect shifts appear which should be removed by a proper rescaling 
of the reduced partition function of CMM. The expression ( |124D reproduces the connection 

The equation ( p.24| ) is based on the following identity for the F-functions 



E c:,^p-ii-r'-'-'-' HE ^ 



s=0 
p-j oo 



a=lka=0 (na - ka)\T (ka + 



X 



r (mh + 



p—i 
p 



II J2 

b=lmt=fb+l (f^b - fb- l)!r (/fe + 



2p-i 
P 



^ mt,=^ ka+i-p+s+l 



1 if fh = Y^ria + n ^- i - p -\- j 
otherwise, 



(125) 



which will be proved in several steps. 
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First, we note that the product of the sums over b in the second hne of (|125|) can be 
easily simplified using the identity 

^ I ^ T{c + b + a-l) 

^=0 7!(a - 7)!r(& - 7)r(c + 7) r(a + l)r(6)r(c + a)r(c + 6 - 1) ' ^ ' 

and the product is equal to 



(127) 



Second, using the identity ( 165D we perform the sum over s in ( |125| ) to obtain 



a=l ka=0 {na - kaY-T (ka + ^) 



, r(Efe. + .-p + i + j + (p-j)^^ 

(EA;a-E/6-p + j + ^ + ^)!r(E/6-ri-p+l + (p-j)^' ■ 



At last, using the identities (|126D j — 1 times and (0), we conclude that ( |125| ) is correct. 
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